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Abstract 

In this paper we introduce a transformation technique, which can on the one hand be used to 
prove existence and uniqueness for a class of SDEs with discontinuous drift coefficient. One the 
other hand we present a numerical method based on transforming the Euler-Maruyama scheme for 
such a class of SDEs. We prove convergence of order 1/2. Einally, we present numerical examples. 
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1 Introduction 

We consider a time-homogeneous stochastic differential equation (SDE) 

dXt = dt + a{Xt)dWt , (1) 

where tr is a Lipschitz, M-valued function and where ^ is an M-valued function that is allowed to have 
discontinuities. 

It is well-known that the problem of existence and uniqueness of a solution to ([T]l is readily settled by 
Picard iteration, if // is Lipschitz, too. In the case where the diffusion coefficient is bounded, Lipschitz, 
and (partly) uniformly elliptic, and the drift coefficient ^ is only bounded and measurable, the pioneering 
work by Zvonkin iflSl and Veretennikov lfT3l [T4ll yields existence and uniqueness of the solution. There 
the result is achieved by applying a transform that removes the drift. This transform can in principle 
be computed by solving a non-degenerate elliptic partial differential equation. In the one-dimensional 
case this reduces to solving an ordinary differential equation. Nevertheless, from the point of view of 
numerical treatment of ([T]) this transformation method is impractical. 

In Leobacher et al. IfTTIl one can find an exisfence and uniqueness resulf for mulfi-dimensional SDEs, 
for fhe case where fhe driff is allowed fo be discontinuous af a hyperplane, or al a hypersurface, buf is 
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well behaved everywhere else. Instead of removing the whole drift by a transformation, only the discon¬ 
tinuity is removed by using a different transformation method. For computing this transform one only 
needs to solve a parametrized family of ODEs which can be done by iterated integration. This transform 
is therefore much less costly from a computational point of view. 

In this paper, we prove another existence and uniqueness result for ([T]l under conditions weaker than 
those in Leobacher et al. m. The transform constructed here is explicit and thus gives rise to a numeri¬ 
cal method that does not require solving any (partial) differential equation. 

In setups with non-globally Lipschitz drift coefficient, various authors have studied convergence of 
numerical schemes. Berkaoui 0 proves strong convergence of the Euler-Maruyama scheme for 
drift. Hutzenthaler et al. ||71 present an explicit numerical method for which they are able to prove strong 
convergence in case of an SDE with non-globally Eipschitz coefficient. Gyongy Q proves almost sure 
convergence of the Euler-Maruyama scheme in the case where the drift satisfies a monotonicity con¬ 
dition. Halidias and Kloeden |ihi] show that the Euler-Maruyama scheme converges strongly in case 
of a discontinuous monotone drift coefficient, e.g., in the case where the drift is a Heaviside function. 
Kohatsu-Higa et al. lITOl show - in case of a discontinuous drift - weak convergence of a method where 
they first regularize the drift and then apply the Euler-Maruyama scheme. Etore and Martinez ||3l 0 
present an exact simulation algorithm for SDEs with a bounded drift coefficient that has a discontinuity 
in one point, but is differentiable everywhere else. 

In contrast to that, we allow the drift to have a finite number of jumps, but require it to be Eipschitz 
otherwise. The transformation method and the scheme based on it are the main contributions of this 
paper. The transformation itself is presented in a constructive way. Moreover, it is chosen such that all 
the involved functions can be computed efficiently. 

The paper is organized as follows. In Section we introduce a transformation method different to 
those from ifTTl [TSll and show that this also leads to an existence and uniqueness result. In Section 
we present the numerical method which is based on the transformation introduced earlier and we prove 
convergence of strong order 1/2. Einally, in Sectionwe present numerical examples. 

2 Existence and uniqueness 

Eet {£, T, {T)t>o, P) be a filtered probability space carrying a standard Brownian motion W = {Wt)t>o- 
Eet /X, fj : M —>• M, be measurable functions. We study the time-homogeneous stochastic differential 
equation (SDE) 


dXt = dt + a{Xt)dWt, Xo = x. ( 2 ) 

The function /x is allowed to be discontinuous. However, the form of the discontinuities is a special 
one: we allow only discontinuities in a finite number m of distinct points < • ■ • < and we assume 
that /X is Eipschitz otherwise. 

Definition 2.1. Eet I C M be an interval. We say a function / : I —M is piecewise Eipschitz if 
there are finitely many points . . . <C ^rn ^ I such that / is Eipschitz on each of the intervals 

(-00,6) n /, (?m,oo) n I and (6,6+i), A; = 1 ,... , 771 . 

We have not assumed anything about the behaviour of / at 6) • • • > However, we have the fol¬ 
lowing elementary lemma: 

Lemma 2.2. For a piecewise Lipschitz function f the one-sided limits always exist. 
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Proof. By symmetry it is enough to show that the left limit exists in every point. For given x G I 
consider a non-decreasing sequence (x„) with Xn x. According to our assumption there is some L 
and some no such that \ f{xn) — f{xk)\ < L\xn — Xk\ for all n,k > no- Now 

OO OO OO 

^ \ f{Xn+l) - f{Xn)\ < ^ L\Xn+l - Xn\ = L {Xn+1 “ Xn) = L{x - Xno) < OO 
n=no n=no n=no 

such that Yl’^=no ifi^n+i) — f{xn)) Converges absolutely and the limit of this series is lim^^oo f{xn) — 
f{Xno)- □ 

This will enable us to transform the SDE into one with Lipschitz coefficients. The transform pre¬ 
sented here is similar to the one in ifTTl but much simpler to compute. 

Assumption 2.3. We assume the following for the coefficients of Q: 

(i) /i is piecewise Lipschitz; 

(ii) cr is globally Lipschitz; 

(iii) There exists c > 0 such that cr^(^i) > c for alH = 1,..., m. 

We want to stress that those assumptions are satisfied by many pracfical examples. A classical one is 
fhe one-dimensional process X satisfying Xq = x and 

dXt = - sign{Xt)dt + dWt ■ 

Our assumpfions on /r enable us fo remove fhe disconfinuify from fhe driff by a suifable fransfor- 
mafion of X. The fransformafion g is chosen such fhaf g" is piecewise linear and is non-zero only on 
environmenfs of fhe disconfinuifies. This yields a piecewise quadratic funclion g', which is consfanfly 
1 excepf on environmenfs of fhe disconfinuifies. Lurfhermore, we can guaranfee boundedness of fhe 
derivatives of g. 

Proposition 2.4. Let with < ■ ■ ■ < and let ai,, am, ..., (3m G M. Let 

0 < >c < 1 be fixed. 

There exist functions g,h '.M —> M such that 

1. g{h{z)) = z for aZZ z G M and h{g{x)) = x for all x G M; 

2. 5,Zi G C^M); 

3. \\g'— iWoo ^ X and \\h'— l\\oo < >c; 

4. supa-gK \g{x) — x\ < ^ max (1, maxi<fc<m(^fc-i-i — ^k)) and g{x) = x for |x| sufficiently large; 

5. g”, h" are piecewise continuous and bounded on M, both one-sided limits of g" at ^k exist, and 
both one-sided limits of h" at g{fik) exist for k = 1,... ,m; 

6- g”{Ck+) = Oik and g”{^k-) = fik, k = 1,... ,m. 

Proof Let ^0 = 6“ ^m-ri = Cm + 1, ao = /3o = am-n = /3m-ri = 0. We construct g” on the 
intervals [Cfc-i, Cfc] for Zc = 1,..., m -|- 1. 
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Choose Cfc E [^k-i,^k-i + + 6)) (cfc will depend on x) and let g" on 5(6-1 + 6)) 

be 


9"ix) 


3a{x—^) 


8a{x—c) 

^ 3(c-?) ’ 

ia{^—c+2{x—c)) 

%=?) 

8«(^—2c+a:) 
3(a;-«) ’ 

lo, 


X E (6 5(^ + 6] 

X E (|(C + c),c] 

X E (c,c+ 

xE(c+£^,c+M] 

xE(c+M,2c-e] 

X E ( 2 c — 6 5(^ + 6)]) 


where we write ^ = 6 -l^ « = o-k-i, and c = Ck for brevity. Note that 


g''{x)dx = 0 and 


ik-l 

and that we choose such that 




max 

ik-l+^k 


'ik-l 


'ik- 


g"{t)dt 


/ g"{t)dt dx = 0 , 
1 J ik-l 


< x/(l + x). 


On [ 2 ( 6-1 + 6 )) 6 ] we define g" analog with g"{^k) = fdk- Further define g"{x) = 0 for x < 6 and 

^ ^ ^m+l- 

Now lef 


pX pt 

g{x) = X + / / g"{s)dsdt. 

Jo Jo 

Then if is easy fo verify fhaf g satisfies ifems 2-6. 

The funcfion g is piecewise cubic wifh positive derivafive and fherefore has a global inverse h fhaf 
can be given explicifly as a piecewise radical funcfion. We have 


\h'iz) - 1 


9'{h{z)) 


l-g\h{z))\ ^ 
\ 9 '{h{z))\ 


□ 

The proof of Proposifion |2.4| is consfrucfive. So now we know how we can consfrucf fhe funcfion for 
fhe numerical approximation. Figure shows fhe functions g' and g” on fhe righf of some 

The funcfion g” is piecewise Lipschifz. The funcfion g', or fhe consfanfs ai, Pi,, am, Pm are 
chosen such fhaf fhe disconfinuifies are removed from fhe driff in a way such fhaf fhe remaining ferm is 
locally Lipschifz, i.e., ak = 2, Pk = 2» and pk = ■ We define fhe 

fransformed SDE by Z = g{X). Then 

dZt = g'dXt + ]^g''d[X\t = dt + ag'dWt. (3) 

Nofe fhaf 0 < x < 1 was arbifrary. Therefore, if x is very close fo 0, fhen g is close fo fhe idenfify. 
However, we wanf fo sfress fhaf fhere is no need fo make x parficularly small. On fhe confrary, in 
general smaller x gives rise lo bigger Lipschifz consfanfs of fhe fransformed coefficienls. x only needs 
lo be smaller fhan 1 fo ensure fhaf g is slricfly increasing. Nofe furlher fhaf around each disconlinuify we 
build Iwo such splines - one on fhe righf hand side, and one on fhe lefl hand side of fhe disconlinuify. 
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Figure 1: The functions g' and g” close to 


This is done to reduce the Lipschitz constant of the coefficients of the transformed equation compared to 
a one-sided compensation of the jumps. 

In Sectionj^we plot the coefficients of the transformed SDE for Example |4.2[ 


Now we are ready to prove the existence and uniqueness result, which is obtained by applying the 
transformation g. 


Theorem 2.5. Under Assumption 2.3 we have that for every x G 
solution Z to the SDE 


there exists a unique global strong 


dZt = p{Zt)dt + d{Zt)dWt , 

Zo = g{x ), 

where fl{z) := p{h{z))g'{h{z)) + ^a^{h{z))g"{h{z)) and cr(z) := a{h{z))g'{h{z)). Furthermore, 
h{Z) is a unique global strong solution to 

Eor the proof we need the following elementary lemma: 

Lemma 2.6. Let f : M —> M be piecewise Lipschitz and continuous. 

Then f is Lipschitz on M. 

Proof. There exist ai < ... < Om such that / is Eipschitz on each of the open intervals (—oo, ai), (ai, 02 ), • ■ •, {am, 00 ). 
Eromthe continuity of / we conclude that / is Eipschitz on each of the closed intervals (—oo,ai], [ 01 , 02 ],..., [om,oo). 

Eet Lq) • • •) Lm denote the respective Eipschitz constants and let L = max(Lo,..., Lm)- Now let 
X, y G M. W.l.o.g. y < X. X and y are in the same interval, then it is obvious that \f{x) — f{y)\ < 
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L\x — y\. Otherwise, we have y < < ... < aj < x with k chosen minimal and j chosen maximal. 


\f{x) - f{y)\ = \ fix) - f{aj) + f{aj) - /(aj-i) + ... + /(a^) - f{y)\ 

< \ f{x) - f{aj)\ + \f{aj) - f{aj-i)\ + ... + |/(ofc) - f{y)\ 

< L\x — Ojl + L\aj — aj-i \ + ... + L\ak — y\ 

= L{x - aj) + L{aj - aj-i) + ... + L{ak - y) 

= L{x -y) = L\x - y\. 


Thus the assertion is proven. 


□ 


Proof of Theorem [23] Define /i(x) := y,{x)g'{x) + \a‘^{x)g" {x) and (t(x) := a{x)g'{x). 

Due to Assumption |2.3| ([II|), a is globally Lipschitz. Furthermore, g' is differentiable with bounded 
derivative, which implies Lipschitz continuity. So a and g' are both Lipschitz and bounded, thus a = ag' 
is Lipschitz. 

To show that ft is Lipschitz, we first note that g is chosen in a way such that fi is continuous. We 
observe that for |x| large it holds that gfx) = 1 and g”{x) = 0. Thus there exists a > max(|^i|, l^ml) 
such that pL is Lipschitz on (—oo, —a), (a, oo). 

Furthermore, f, is Lipschitz on the intervals (—a, ^i), (Ci, ? 2 )j • • •, Cm), (Cm, a) as a sum of 

products of bounded Lipschitz functions. That means that fl is piecewise Lipschitz. 


Hence f, is Lipschitz by Lemma 2.6 


Since g' is bounded away from 0, h' is bounded and thus h is Lipschitz. Thus, also jl, a are Lipschitz. 


From the Lipschitz continuity of f, and a we get existence and uniqueness of a global strong solution 
to Q from ifT^ Theorem 3.1]. Furthermore, note that Ito’s formula holds for /i by lUl Problem 7.3]. By 
applying Ito’s formula to h we get that there exists a unique global strong solution to □ 


3 Numerical scheme 

The numerical scheme is based on the transformation g introduced in Section In fact, we transform 
the initial value by applying g and then solve the transformed SDE by applying the Euler-Maruyama 
method. This converges strongly with order 1/2 to the solution of the transformed SDE. Then we apply 
h to get the solution of our original SDE. 

Eor any f > 0 and <5 > 0 we denote the n-step Euler-Maruyama scheme with step-size <5 recursively 
by 


t, 5) ■.= z + p.{z)6 -h a{z) {Wt+5 - Wt) 
t, 6) := (j)'^{(f)^{z, t,5),t + 6,S). 

The scheme for solving @ looks as follows: 

^{x,t,T,n) := h{f'^{g{Xt),t,T/n)). (4) 

Now we are ready to prove convergence of scheme 0 - 

Theorem 3.1. Let Assumptions ^f3\ hold. 

Then scheme 0 converges with order 7 = Xj^lto the solution of 0 , i.e., 
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Proof. 


E {\\Xt - Hx,0,T,n)ff" = E iWHgiXT)) - hir{g{x),0,T/n))\Y" 

<LhE{\\ZT-rigix),0,T/n)ff^ , 

where Zt is the exact solution to the transformed SDE. From |l9j Theorem 10.2.2] and as /2 and a are 
Lipschitz we have 

E{\\ZT-r{z,0,6)f) <Cd. 

Altogether this yields 

E{\\XT-Hx,0,T,n)\\Y^ <C6^/\ 
where C = and 6 = T/n. 

□ 


4 Examples 


In this section we present numerical examples, where we show the exact choice of the transformation and 
the transformed parameters. Furthermore, we investigate the convergence of our method and compare it 
to the Euler-Maruyama method. 

Example 4.1. First, we consider the process X mentioned in the beginning and satisfying the equation 

dXt = - sign{Xt)dt + dWt ■ 

That is, g,{x) = — sign(x) and a{x) = 1. We choose different values of x = 1/16,1/64,1/256. 
Figure|^shows the estimated £^-error between two consecutive discretizations of the transformed Euler- 
Maruyama method (EMT 1/x) in comparison to crude Euler-Maruyama (EM). This means we calculate 



x!^'^ — x\^ ] and plot it over log(i(*^\ where xl^^ is the numerical approximation 


with stepsize 5 = 5^^'^ and E is an estimator of the mean value using 1024 paths. 

We see that for this parameter choice the crude Euler-Maruyama method is even better than the 
transformed method. Furthermore, we see that the estimated -error for the crude Euler-Maruyama 
method seems to be of order 1. The reason for this is probably that in case of Eipschitz coefficients 
Euler-Maruyama converges with strong order 1, if the diffusion parameter is constant. However, we 
cannot expect strong order 1 for the transformed method, since a is not constant, even if a is constant. 

In the next example g has more than one discontinuity. In addition to the estimated £^-error between 
two consecutive discretizations we will show how the transformation looks like. 

Example 4.2. Now, we consider the SDE 

dXt = g{Xt)dt + a{Xt)dWt . 

The coefficients are chosen as follows: 

X — 2 , X < —\ 


g{x) = < 


2 , 


-1 < X < -0.5 


1 — , —0.5 < X < 0 


x^ , 0 < X < 1 

^ —X — 1 , X > 1 , 


7 






Figure 2: The estimated £^-error for different values of x for Example |4.l[ 



Figure 3: The functions jj. and a. 


Figure shows the parameters /x, a. Figure shows the derivatives of the transformation g, and 
Figure]^ shows the parameters fl,a of the transformed SDE. 

Again, we choose different values of x = 1/16,1/64,1/256. Figure shows the estimated C^- 
error between two consecutive discretizations of the transformed Euler-Maruyama method (EMT 1/x) 
in comparison to crude Euler-Maruyama (EM). 

In this example we reach the calculated order 1/2. Whether the transformed method is better than 
Euler-Maruyama depends on the choice of x. 

As the last example we solve an SDE appearing in insurance mathematics. 

Example 4.3 (Threshold dividend strategy). In |[I] the authors study the dividend maximization problem 
from risk theory in a diffusion model. They find that the optimal dividend policy is of threshold type 
with constant threshold level b. This means that dividends should be paid at the maximum rate whenever 
the surplus process X of the insurance company exceeds b. Otherwise, no dividends should be paid. Our 
numerical scheme enables us to simulate the optimally controlled surplus process from |T) : 

dXt = [9 - K f{x>b}) dt + adWt , 

where 9 is the drift of the uncontrolled process, and K is the maximum dividend rate. We choose 
9 = 1, K = 1.8, CT = 1, and b « 0.895635 is the optimal threshold level, which can be calculated as in 

m. 
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Figure 5: The transformed parameters fl (blue) and a (red). 



Figure 6: The estimated £^-error for different values of x for Example |4.2[ 
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Figure 7: The estimated £^-error for different values of x for Example|4.3[ 


Figure 1^ shows the estimated -error between two consecutive discretizations of the transformed 
Euler-Maruyama method (EMT 1/x) in comparison to crude Euler-Maruyama (EM) for x = 1/16,1/64,1/256. 

Threshold type control strategies appear frequently when solving stochastic optimization problems 
in various fields of applied mathematics. Therefore, our scheme potentially serves for a wide range of 
applications. 
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